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The IKKT or IIB matrix model has been proposed as a non-perturbative definition of type IIB superstring 
theories. It has the attractive feature that space-time appears dynamically. It is possible that lower dimensional 
universes dominate the theory, therefore providing a dynamical solution to the reduction of space-time dimen- 
sionality. We summarize recent works that show the central role of the phase of the fermion determinant in the 
possible realization of such a scenario. 



One of the most important problems in string 
theory is the emergence of four dimensional 
space-time in low energy physics. Compactifi- 
cation schemes have been the most popular ap- 
proach of the issue, but the price to pay is loss 
of predictability. One possible scenario is that 
4d space-time appears dynamically from an orig- 
inally higher dimensional theory, usually 10-D or 
11 D. A promising model that might realize this 
scenario is the IKKT matrix model which is 
supposed to define type IIB superstrings in the 
large N limit non-perturbatively. One can view 
this model as a "lattice string theory" and study 
it using standard Monte Carlo techniques. For- 
mally it is the zero- volume limit of-D = 10, J\f = 1 
super Yang-Mills theory. The partition function 
of the model (and its generalizations to D = 4, 6) 
can be written as 

Z = J dAc- Sb Z t [A] , (1) 

where (fi = 1, •••,£)) are D bosonic N x 
N traceless hermitian matrices, and Sb = 
— jp-Tr {[Ap, A u ] 2 ) is the bosonic part of the ac- 
tion. The factor Zf[A] represents the quantity 
obtained by integration over the fermionic ma- 
trices, and its explicit form is given for example 
in Ref . . Space-time appears dynamically (the 
eigenvalues of the bosonic matrices A^ represent 
space-time points) and if the dominant configura- 



tions have d extended dimensions and D — d small 
ones then we obtain essentially a d dimensional 
space-time dynamically. This requires in partic- 
ular the spontaneous symmetry breaking (SSB) 
of the manifest SO(D) invariance of the model. 
We stress here that the parameter g appearing 
in eq. (|l|) is not a coupling constant but a scale 
which can be absorbed by a redefinition of the 
fields. One hopes to arrive at d = 4 in the 1QD 
IKKT model. 

Such a scenario has not been verified. There 
are severe technical problems in studying the 101? 
IKKT model by computer simulations, related to 
the complex action problem. The study of sim- 
pler models has shed some light on the possible 
mechanisms that could realize the above men- 
tioned scenario. We have learned that SSB is 
not realized in the absence of the rapidly oscillat- 
ing phase r of the (complex) fermionic partition 
function Zf [A] . In Ref. || it has been shown that 
when fermions are absent ( "bosonic model" ) SSB 
does not occur. The study of a low energy ap- 
proximation of the 101? and 6D models, where T 
is set to zero by hand, did not reveal any indica- 
tion for SSB {§. When the effect of T is infinitely 
enhanced in the 10-D case by appropriately de- 
forming the original model, it turned out that 
space time is 2 < d < 8, i.e. the SSB scenario 
is realized |5) . Recently it has been shown in ex- 
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actly solvable matrix models that SSB of SO(D) 
occurs precisely due to the phase of the fermion 
determinant |J. 

The purpose of this work is twofold: First to 
show that in the AD model SSB does not occur 
@ . The AD model has the property that Zi [A] is 
real positive, T = 0, and it can be studied using 
ordinary Monte Carlo methods. Therefore the 
result supports the importance of the phase in 
the realization of SSB. The authors of Ref. Q 
have raised the possibility that one dimensional 
structures dominate eq. (|l|). We propose differ- 
ent, more physical probes of space-time dimen- 
sionality that need to be used in D = A. Due to 
the power tail of eigenvalue distributions JjJ, the 
observable = -^Tr (A^A^) used as a probe 
in Ref. g is ill defined (note that (Tr A 2 ) is di- 
vergent). This does not happen in higher dimen- 
sions where one expects all probes discussed here 
to give identical answers. Of course in any dimen- 
sion (Tr A m ) diverges for large enough m, there- 
fore it is important to understand and resolve this 
puzzle in a convincing way. 

Second we show that it is possible to attack the 
problem of SSB in higher dimensions despite the 
complex action problem. We propose a method 
that allows us to compute the T^ v eigenvalue dis- 
tribution in 61? for large N. The method does 
not suffer from the complex action problem and 
it might be applicable to other models as well. 
The resulting distribution is qualitatively differ- 
ent from the case T = 0, due to the presence of 
the phase. A double peak structure appears for 
large enough N, raising the possibility that small 
eigenvalues dominate for some dimensions while 
large ones dominate for others, which would mean 
that SSB occurs. A more involved study, possible 
however, is needed in order to resolve the issue. 

A natural probe of space-time dimensional- 
ity is its "moment of inertia" tensor = 
jjTr (A^A V ). Its D eigenvalues Ai > A 2 > . . . > 
Xd > represent the principal moments of in- 
ertia. In AD, however, its largest eigenvalue is 
known to diverge 0, in agreement with the di- 
vergence of (Tr A 2 ) ||. In this case it is possible 



to modify its definition to 

y(new) _ / 2( x ifi ~ Xj[i){xii> ~ %jv) \ 

~\f^ N(N-iW( Xi - Xj y I ' 

are the space-time points defined to be the di- 
agonal elements of the bosonic matrices A^ when 
they are transformed into a form as close to si- 
multaneously diagonalized as possible . 
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Fig. 1: The eigenvalues x\ ncw) /(^/gN 1 / 4 ). 
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Fig. 2: The four Polyakov lines (P tl (p)) (/J, = 
l,...,4j atp 2 = 0.2. 

Since the distance between space-time points 
behaves asymptotically as p(r) ~ r~ 3 , this quan- 
tity is well defined. In Ref. fjj we found that 
for increasing N, the eigenvalues / \( new ) , . . X^ ew ^ 
show no trend for SSB. This can be seen in Fig. 
1 where all eigenvalues seem to approach the 
same value with increasing N. The same re- 
sult is obtained from the study of Polyakov lines 
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P(p) — -kTi exp(ip A1 A M ). They are essentially 
the Fourier transform of p(r) and their character- 
istic fall off is a measure of the extent of space- 
time Q . Since Wilson loops in general have been 
interpreted as string creation operators they cor- 
respond to observables in string theory. Defining 
the Polyakov lines on field configurations that 
diagonalize as P^{p) — jj Tr exp(ipA^) we 
can study if SSB occurs. In Fig. 2 we see that 
all Pfiip) seem to converge to the same univer- 
sal function, so no trend for SSB is observed 0. 
This behavior is typical for all values of p 2 that 
we measured. 

For the 6D partition function eq. ([[]) T ^ and 
the model suffers from the complex action prob- 
lem. The partition function can be written as 
Z = JdA e- So e lT , where Z { [A] = exp(r R + ir), 
r R , T G K and S = S b - T R . The rapidly fluc- 
tuating e lF term makes the calculation of expec- 
tation values exponentially hard with increasing 
system size. Furthermore, the model defined by 
the partition function Z a = J dA e~ 5 ° used in the 
Monte Carlo simulation visits very rarely the part 
of configuration space that dominates in the full 
model (overlap problem). It is possible though 
to overcome these technical obstacles. One can 



study the distribution of Xi 



(\i)o 



defined by 



Pi(x) = (S(x — Xi)), where (. . .) are expectation 
values with respect to Zq. In Rcf. |I(J we observed 



(2) 



where pf (x) is the Xi distribution in the 
Z ensemble, = log (cosr) ^. =a . 

and C — (cosr) is a normalization con- 



stant. 



are expectation values in the 

Z Q ^. =x = f dAe~ s " 5(x — Xi) ensemble. It turns 
out that the function <5>i(x) converges to an N- 
independent scaling function for quite small N 
(< 20). Combining this observation with the fac- 
torization property (Q), we see that pi(x) can be 
computed for much larger N than (x) since the 
remaining terms in eq. (||) do not suffer from the 
complex action problem. 

In Ref. 1 10 1 we simulated the low energy ver- 
sion of the 6D IKKT model for N up to 128. 
We computed the scaling function $i(x), i = 4, 5 



using N < 20. We found that pi{x), i = 4,5 
have two peaks at x — x< and x — x> (a;< < 1 
and x> > 1), which is qualitatively different from 

Pi°\x), with peaks around 1. It is important to 
determine the peak which dominates the large- 
N limit. In Ref. (h]] we found that the data are 
consistent with a scenario where a;< dominates 
p$(x) whereas x > dominates pi{x). This behavior 
can be expected since the qualitative difference 
in the behaviors of ps,(x) and pa{x) is mainly due 
to p^\x) and pf\x). In the branched polymer 
description of the low energy effective theory of 
the IIB matrix model |l| , pf^ (x) is expected to be 
much more suppressed in the small x regime than 
p^\x). The method described above is quite 
general and can be applied to other physical sys- 
tems that suffer from the complex action problem. 
Such a system is finite density QCD, and we are 
currently investigating the related random matrix 
theory with very encouraging results jOj. 
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